AIP/123-QED 



Current-induced phase transition in ballistic Ni nanocontacts 

R.G. GatiyatovH V.N. Lisin, and A. A. Bukharaev 

Zavoisky Physical- Technical Institute RAS, ^20029 Kazan, 
Russia 

(Dated: 4 March 2010) 

Local phase transition from ferromagnetic to paramagnetic state in the region of the 
ballistic Ni nanocontacts (NCs) has been experimentally observed. We found that 
contact size reduction leads to an increase in the bias voltage at which the local phase 
transition occurs. Presented theoretical interpretation of this phenomena takes into 
the account the specificity of the local heating of the ballistic NC and describes the 
electron's energy relaxation dependences on the applied voltage. The experimental 
data are in good qualitative and quantitative agreement with the theory proposed. 
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For more than a decade a point contact between two metal conductors with characteristic 
micro- and nanoscale dimensions is one of the interesting objects studied theoretically and 
experimentally^. Recently, the investigations were mainly concentrated on the atomic scale 
contacts due to the rich variety of the quantum size effects observed in such systems^-. 
However, the nanocontacts (NCs) with the transverse dimensions of the order of an electron 
Fermi wavelength are not suitable for practical usage because of their rapid destruction. 
Therefore, ballistic NCs with the diameter d larger than Xp, on the one hand, but smaller 
or comparable to the transport electron mean free path l tr , on the other hand, are more 
promising for applications. The investigation of the specificity of heating of the magnetic 
ballistic NCs and the adjacent regions by electrons current is one of the topics that is 
particularly due to the spin transfer torque effect in such structures 6 where temperature can 
play a crucial role^. Moreover, the heating of the contact region to the Curie temperature 
should lead to a phase transition from ferromagnetic to paramagnetic state. Earlier, such 
phase transition was achieved in a number of ferromagnetic microcontacts (ltr « d) by 
current-induced heating at liquid helium temperatures 8 . However, up to now the phase 
transition of the ferromagnetic ballistic NCs was not investigated. In this case (d ~ lt r ), 
it is known that the voltage drops in the region of the order of the contact diameter and 
the electrons release their excess energy beyond the the region of the potential drop^. It 
determines the specificity of the thermal heating of the contact and the adjacent regions 
that is of interest of this research. 

Ni NCs has been formed between two microwires fixed on a substrate using electrochem- 
ical method 10 . NCs were fabricated in the nickel sulphate solution 0.25 M MSO4 + 0.5 M 
H3BO3 (working voltage 1-1.4 V). The conductance and the current-voltage (I-V) charac- 
teristics were recorded with two digital multimeters Agilent 34410A using 4-probe method. 
A single 70 Hz triangular voltage pulse was applied to the circuit to obtain current-voltage 
curves. I-V curves were recorded at room temperature in the bath with bidistilled water, 
which has the conductivity much smaller than a NC conductivity. The zero-bias resistance 
of the fabricated Ni NCs was in the range of 30-400 Ohm. 

A typical I-V curve and the dependences of the dR/dll spectrum and the resistance of 
the Ni NC versus applied voltage are shown in Fig. [U Forward and backward branches of 
the current-voltage curve coincide. Therefore, we conclude that the current does not affect 
the contact, and softening effect 11 or influence of the solution^ are absent or negligible in 
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FIG. 1. I-V curve (top) and typical dependence of the resistance and dR/dU spectrum of the Ni 
NC versus applied voltage (bottom). 

our case. 

I-V curves are nonlinear and the resistance increases with the increase in the applied volt- 
age. Moreover, the dR/dU spectra have a maximum. The position of the maximum depends 
on the resistance (size) of the Ni NC (Fig. [2]). It is worth noting that such experiments were 
also carried out for Cu NCs. In case of Cu the dR/dU spectra do not contain any features 
in the same region of the applied voltages. 

The observed feature on the dRj dU spectra we attribute to the local phase transition from 
ferromagnetic to paramagnetic state in the vicinity of the Ni NC region. Above the Curie 
temperature (Tc — 631 K for Ni) the electron-magnon scattering saturates and becomes 
temperature- independent 12 . So, dR/dU is determined only by electron-phonon scattering. 
It should lead to the maximum in the dR/dU spectrum versus U that was observed experi- 
mentally (Fig. [p. 

From the solution of the steady-state heat conduction equation one can find the contact 
resistance Rc at critical voltage Uc at which the local phase transition occurs . Let us 
make several assumptions. First of all, let us assume that the form of the NC matches with 
an orifice that is the diameter of the NC d is much larger than its length. Secondly, let's 
assume that the relaxation of the electron's excess energy occurs in a sphere with the radius 
b around the contact. Assuming the density of the released thermal power constant, the 
problem becomes spherically symmetrical. 

In the steady state, the thermal flux through the sphere with the radius r is equal to the 
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FIG. 2. Dependence of the bias voltage Uc at which the local phase transition occurs on the 
resistance Rc of the Ni NCs at this voltage (see Fig. [T|) . Each dot was obtained from the individual 
dR/dU spectrum of the Ni NC with different zero-bias resistance. Solid line is the result of the 
fitting by Eq.flTJ). 

released power in this sphere: 



where q = — AVT is the thermal flux, A is the volume thermal conductivity, U is the voltage 
applied to the NC, R is the NC resistance, S r is the surface of the sphere with radius r. By 
solving the system of the equations (OQ) one can obtain the relation between the temperature 
Tf, on the sphere with the radius b and the temperature T nc in the center of the NC (r = 0): 



where T is the the temperature of the conductors far away from the NC. 

The b value is determined 9 by the dimension of the region where the external potential 
drops (of the order of the NC diameter- 1 ^) and the electron diffusion length over the 
energy relaxation time t e and can be estimated as: 




(1) 
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r, 





b = d + A E , 



(3) 



where 



A E = y/6D tr (T)r E (T, U) = ^2l tr (T)l E {T ) U). 



(4) 
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Here D tr = \vplt r is the transport diffusion coefficient of the electrons, vp is the Fermi 
velocity, l tr = r tr vp, lp = T E vp, r tr is the transport relaxation time. 

Now we take a look at the mechanisms of the momentum and energy-loss mechanisms 
of the hot electrons in more detail. In ferromagnetic metals, the excess energy E of the 
electrons can be transferred either to the lattice via the phonon emission or to the magnetic 
subsystem via the magnon emission. If the excess energy of the electron E strongly exceeds 
the energy of the phonons with the Debye frequency (Ep> = 35meV for Ni^), then at each 
act of the phonon emission an electron radiates a phonon with the Debye frequency and 
undergoes large- angle scattering^. Therefore, only one act of the scattering is needed for the 
electron in order to forget the direction of its motion. The mean time between such emissions 
r ep h determines the electron-phonon contribution to the transport relaxation time Tt r . For 
the time t the electron emits t/t ep h phonons with the Debye frequency and looses energy 
Ep,t/T ep f l . In the case of the electron-magnon scattering, the mean time between emissions 
of the magnons r em plays the role of T ep h and the exchange energy E ex (E ex = 54.5 meV for 
Ni^) plays the role of Ep. Then the effective time Tp needed the electron to entirely lose 
the excess energy E is proportional to this energy: 

E (5) 

Ed/ T eph + E ex / T em 

From Eq. (jSj) one can see that for E » Ep>, E ex the energy relaxation time Tp becomes much 
larger than the transport time r tr (neglecting electron-impurity scattering r^ 1 = T~ ph + 
)• 

Let us now consider the case d < lt T - In this regime, the electron ballistically passes 
the region where most of the applied voltage drops and it gains the energy E = \eU\ (e is 
the electron charge). By substituting Eq. (JSJ) into Eq. (J3J) and taking into the account that 
d << Ap, the expression for b is: 



b(f,U) ~ \l gsgM! ^ K y^, ( 6) 

Ep/hphiT) + E ex /l em (T) 

where T is the temperature averaged over a sphere with radius h (the temperature is not 
constant inside the sphere), l ep h = r ep hVF, l em = r era vp. Therefore, according to Eq. (JHJ), 
the size of the heating region b depends on the applied external potential difference. This 
effect separates the ballistic NCs from the microcontacts and determines the specificity of 
the current-induced heating. 



By substituting Eq. ([6]) into Eq. ([2]), one can obtain the dependence of the voltage Uc 
needed to heat the NC to the critical temperature Tc on the NC resistance Rc'- 

U c = AR C , (7) 

where A is constant: 
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(8) 

Let us now demonstrate that Eq. (j7j) can be applied to the experimental dependence 
Uc{Rc) i- e -> the condition d < l tr (Tc) holds. The diameter of the fabricated Ni NC at 
To = 300 K can be found using the value of its zero-bias resistance from Wexler formula^. 
To this end, one needs to know lt r {To) and the value of the electrical resistivity p(Tq) = 
7.2 • 10~ 8 Ohm-rn^ for Ni. We estimated the values of the mean free paths as l(T) = 
5 • 10~ 16 /p(T). This gives k r {To) = 7nm and l tr (Tc) = 2nm. According to Wexler formula, 
the NCs with the zero-bias resistance of 30, 100 and 400 Ohm have the diameters of 6.4, 3.2 
and 1.5 nm, respectively. Therefore, the ballistic transport of electrons (d < k r {Tc)) should 
exist for the contacts with R > 100 Ohm even at the Curie temperature. 

Another important condition of the theory proposed is the constancy of the phase tran- 
sition temperature for the contacts of different size. We suppose that the surface and size 
effects are negligible and can not strongly affect the value of the Curie temperature of the 
contact region. This means the phase transition occurs at the same temperature for the 
nano contacts of different size. 

A solid curve in Fig. [2] represents the result of the fitting of the experimental data using 
the function Uc — AR C by root-mean-square method. 

The exponential value a = 0.67 ± 0.02, which is in good agreement with the value from 
Eq. (j7|), and the coefficient A = (0.95 ± 0.12) • 10 -2 V/Ohmi were extracted using fitting. 
One can see that the exponential law Eq. ([7]) describes well the experimental dependence 
Uc(Rc) for \eUc\ > Er>,E ex in wide range of the NC resistances. The estimate of the A 
using Eq. (JSJ) gives the value of 2.2 • 10 -2 V/Ohm§, which is of the same order of magnitude 
as the experimental one. Here we used the parameter values T = (Tc + T),)/2 = 570 K 
(T b = 508 K was found from Eq. (j2]) for T nc = T c and the known dependence A(T)^), 
ltr(T) = 2.6 nm, l ep h{T) = 3.9 nm and l em (T) = 8.2 nm. The values of the mean free paths 
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were estimated as l(T) = 5-10 16 /p(T) where the p ep h, p em and p tr = p ep h + Pem values for 
Ni were taken frorn^. 

For the Ni NC with the smallest zero-bias resistance of 30 Ohm (Rc = 40 Ohm) the 
electron transport is diffusive because d > l tr (Tc). The maximum on the dR/dU spectrum 
is observed at 0.12 ± 0.02 V. That is in good agreement with the voltage of 0.19 V needed 
to heat the diffusive Ni microcontact to the critical temperature^. Here one has to take into 
the account that in^ the heating was carried out at the liquid helium temperature. 

The estimates are in good qualitative and quantitative agreement with the experimental 
results. It confirms that the observed feature on the dR/dU spectra is due to the phase 
transition from ferromagnetic to paramagnetic state in the region of the ballistic contact. 

The main specificity of the current-induced heating of the ballistic NCs is in the voltage 
dependence of the heated region size for the high applied voltages (\eU\ » Ed,E sx ). It 
results in the increase in the bias voltage needed to heat NC to the critical temperature 
(Fig. [2]) with the reduce of the Ni contact's size (the increase in the NC resistance). While 
in microcontacts it does not depend on the contact size^. 

It must be noted that I-V curves of the ballistic Ni NCs were studied in the past^ for 
U « Uc, where only a negligible deviation from the linearity was found. The authors 
concluded*^ that I-V curves of clean metallic ballistic Ni NCs must be linear and the origin 
of the nonlinearity is due to the contamination of the contact^. In our case (fig. 2) the 
nonlinearity of the I-V curves is also negligible at small applied voltages as in 16 . While the 
nonlinearity becomes noticeable at high voltage (U ~ Uc)- Uc increases with the reduce of 
the NC size. In case of contaminated NCs the resistance drops with an increase of the applied 
voltag ell that is in contrast to our results. We think the influence of the contaminations 
can be eliminated in our case. 

In conclusion, we have shown that I-V curves of the ballistic metallic Ni NCs (d < k r ) 
are nonlinear in the region of high applied voltages (U ~ Uc) as a result of Joule heating of 
the contact region. Such heating results in the phase transition in the contact region if the 
critical temperature reached. 
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